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Abstract 

This Letter describes the first application of entanglement-induced virtual-state 
spectroscopy to a molecular system. Non-classical, non-monotonic behavior in a two- 
photon absorption cross section of the OH A-X system, induced by an entangled biphoton 
state is theoretically demonstrated. A Fourier transform analysis of the biphoton cross 
section permits access to the energy eigenvalues of intermediate rovibronic states with a 
fixed excitation photon energy. The dependence of the Fourier spectrum on the tuning 
range of the entanglement time T e , and the relative path delay r e is discussed. Our 
analysis reveals that the implementation of molecular virtual-state spectroscopy for the 
OH A-X system requires the tuning of r e over a pico-second range with femto-second 
resolution. 


The unique characteristics of entanglement-induced multi-photon spectroscopy 


which arise from the fourth-order quantum interference of an entangled photon with a 


medium (atom or molecule) have been recently discovered [1-4]. In entanglement- 
induced two-photon laser induced fluorescence, the fluorescence emission from the 
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medium being probed automatically acts as both a detector and temporal coincidence gate 
for the two-photon absorption process [5]. Thus, the entangled two-photon experiment 
would indeed display quantum interference effects in the probed medium without 
employing the requisite two separate detectors and coincidence circuit. The entangled 
two-photon absorption spectrum of atomic hydrogen was explicitly demonstrated 
numerically, and the resulting non-classical behavior such as a linear dependence of the 
two-photon absorption rate on the photon-flux density, and the non-monotonic behavior 
of the absorption cross section versus entanglement time were theoretically described by 
Fei et al [1]. Furthermore, a remarkable feature of this emerging spectroscopic technique 
was illustrated by Saleh et al. [2]: the ability to access intermediate states (or virtual- 
states) of the two-photon transition via a “transition-path-scan” without changing the 
wavelength of the excitation light source. In this Letter we extend this virtual-state 
spectroscopy to a molecular system. 

The hydroxyl radical (OH) was chosen for our demonstration of entangled two- 
photon molecular spectroscopy because OH is one of the simplest molecules, yet is a 
very important species in the chemistry of many processes. The OH molecule is a key 
intermediate in many chemical reaction processes involving oxygen or water such as 
combustion or atmospheric chemistry. The spectroscopic measurement of OH provides a 
way to probe flame chemistry or to trace the reaction pathways in atmospheric chemistry. 
The use of OH in the present study has a further advantage since the molecular structure 
is well known, and its spectroscopy over a wide range of spectral regions has been 
extensively investigated [6,7]. The practical convenience of generating OH through a 
microwave discharge of water, or as provided from combustion products in flames is also 
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an important factor in an effort to implement a proof-of-concept experiment of the 
present technique in future. 

Because of its significance, optical diagnostic techniques that probe OH have been 
thoroughly studied including laser absorption spectroscopy of the ground X state or laser- 
induced fluorescence of the A-X electronic system [6-8]. In particular, two-photon- 
excited fluorescence measurements have been successfully applied to the electronic 
system of OH [9,10] despite its very small cross section. One of the major problems in 
the two-photon spectroscopy of OH, is that it requires a very large excitation photon flux 
to gain sufficient fluorescence signals since the classical two-photon excitation occurs 
with randomly arriving photons (as evident by its quadratic absorption rate). Thus it 
suffers from a higher probability of photochemical reactions in the probe volume formed 
by the intense laser pulses [10,11]. This problem is significant not only for the current 
system but also for other kinds of classical two-photon spectroscopy, especially that as 
applied to the imaging of delicate biological tissues [4], However, by using entangled 
biphoton pairs, the two-photon absorption rate could in theory, become comparable with 
that produced by classical light source but at a far lower excitation photon flux [1,4]. 

This arises mainly because of the twin photon nature of the entanglement process that 
produces a highly time-correlated photon pair that has a high degree of fourth-order 
coherence [12,13]. Thus the risk of photochemical effects can be greatly reduced. Hence, 
entangled biphoton spectroscopy can increase the two-photon absorption efficiency of 
OH, and furthermore allow us to access intermediate state information in an 
unprecedented way. 
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Absorption of an entangled photon pair . — A photon pair produced by the process of 
spontaneous parametric down-conversion (SPDC) within a nonlinear crystal of finite 
thickness, with second-order susceptibility, has the form of the twin state that is 
expressed as an expansion in the wave-vector space [12-14]: 


where k indicates the wave vector and the subscripts 1, and 2 indicate signal, and idler, 
respectively. v F(k,,k 2 ) is a nonseparable wave function in k space that is a function of 

the phase matching condition, i.e. kj +k 2 = k p where p indicates pump. The entangled 

photons are generated via the SPDC process within an entanglement area A e and within 
an entanglement time T e [1,2], which depend on the property of the crystal and the pump. 

We now consider a two-photon transition, where a quantum entangled photon pair 
interacts with a molecule in an initial state | y/j, with an energy eigenvalue e,. The photon 
pair excites the molecule to a final state | y/j) with energy eigenvalue ^through an 
intermediate state | y/j) with energy eigenvalue £j (assuming that an intermediate state 
linewidth can be neglected [1]). Further assuming a monochromatic pump source with 
the frequency co p = £f- based on second-order time-dependent perturbation theory, the 
entangled two-photon (or biphoton) absorption cross section may be given by [1,3] 


where h is the Plank constant, s 0 is the vacuum permittivity, c is the speed of light, and 
of + of = o) p (phase matching condition). Here, D } represents the transition matrix 



0 ) 
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elements for the two-photon excitation process in the molecule, and is given by Dj= {y/f 
\d\ y/j) < y/j\d\ y/t), with rovibrational or electronic transition moment d. The energy 
mismatch is identified as 


where the superscript k = 1 or 2 refers to the signal and idler photons [see Fig. 1], The 5- 
function for the frequency terms ensures that energy conservation is satisfied. When the 
final state of the molecule is assumed to have Lorentzian broadened energy distribution 
with radiative lifetime (or natural lifetime), the cross section cr e may be multiplied by the 
final state linewidth function and integrated over the final-state energy [2], In this case, 
the ^-function removes the integral, and is simplified to become the linewidth function 
through the so-called sifting property. Note that the entangled two-photon absorption 
occurs via a ft/we-entangled quantum state (twin state) given by Eq. (2) which contains T e , 
a variable with time units. This entanglement time dependency does not appear in an 
equation describing a classical two-photon absorption process. 

In order to realize virtual-state spectroscopy, another optical parameter, the relative- 
path-delay T e (| r e | < T e ), has to be externally introduced between the two twin photons 
generated from the SPDC process. Based on Eq. (2) and the concept of integrating the 
cross section over a range of T e , Saleh et al [2] derived the weighted-and-averaged 
entangled-photon cross section for the degenerate pump case as 



(3) 


— 

<7e^— £ -S(£ f -£ i -G) )s(T e ) , 
164 , 
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where 
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s(?') = +^K| 2 [2 + exp(— i'2A jT e )+ exp(i'2A y r e )] . (5) 

j 

Here, the coefficient 4/ = Djl A j. In Eq. (4) the averaging (integration) range of T e has to 
be much greater than the inverse of the smallest frequency difference between 
neighboring intermediate states, i.e. 1/| £) - Sk |. Units here are defined as atomic units and 
such that h = c = 1 . As seen in Eq. (5), the relative phase of the entangled photon pair can 
be tuned by adjusting r e , and this relative-phase tuning varies the absorption amplitude. 
The Fourier transform of a e over a range of r e can then, in principle, reveal the energy 
mismatch terms Ay (in frequency space), thus the intermediate state energy as well as the 
transition probability term Aj (magnitude of the Fourier transform) of the molecule. The 
practical implementation of virtual-state spectroscopy is briefly described by Saleh et al 
[ 2 ]. 

Two-photon electronic transition in A-X system of OH. — To theoretically 
demonstrate entangled two-photon spectroscopy of OH, the i?i(5) (in standard ^AJf/ 
notation, with redundant labels suppressed [8,9]) rotational line of the two-photon A 2 Z + - 
X FI, (0,0) transition is studied. More specifically, this is the rovibronic transition from 
the initial state X 2 Tly 2 (v" = 0, N"= 5) to the final state A 2 Z + (v' = 0, H- 6). We 
throughout use Hund’s case (b) rotational quantum number indication (N) while the 
molecular structure of OH starts as case (a) at low J and switching shortly thereafter to 
Hund’s case (b) due to the spin-orbit uncoupling. Here, the vibrational quantum number 
is denoted by v. The two-photon i?i(5) line has been experimentally confirmed via 
classical two-photon fluorescence technique by Crosley and Smith [9]. 
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The two-photon excitation of OH described above essentially involves pairs of 
transitions that consist of combinations of one rovibrational (infrared) transition and one 
rovibronic (UV) transition. Such transitions can be: 1) an X-X infrared transition 
followed by an A-X UV transition, or 2) an A-X UV transition followed by an A-A 
infrared transition. The selection rule for a two-photon electronic transition is given by 
AN = 0, ±1, ±2, ±3 (A J = 0, ±1, ±2), which dictates the relationship between an initial 
state and a final state. The R\{5) line (AJ = AN = 1) obeys this rule. Since a two-photon 
electronic transition occurs via two single-photon excitation processes, the single-photon 
selection rules also limit each step that is either from an initial to an intermediate state, or 
from an intermediate to a final state. The single-photon selection rule for an X-X 
transition is AN = 0, ±1, ±2 (AJ= 0, ±1). Thus only P, Q, and R branches are allowed in 
2 n 3 / 2 - 2 n 3 /2 transition and Q P, R Q, and S R branches in 2 IIi/ 2 - 2 n 3 /2 transition (hereafter the 
branches are given in ^ AJ notation for rovibrational transition, with redundant labels 
suppressed). As for an A-A transition, the single-photon selection rule is AJ- AN = ±1 
since the total electronic orbital angular-momentum quantum number A = 0 for £ state. 
Thus only p P or r R branches are allowed for this transition. Here satellite branches (AN 
= 0) in an A-A transition is neglected because their intensities decrease rapidly (as an 
exponential function) with N and becomes very weak over N = 5. The single-photon 
selection rule for a A-X electronic transition is AN— 0, ±1, ±2 ( AJ= 0, ±1) and there are 
12 selection-rule allowed branches as described in conventional single-photon induced 
electronic transitions [15]. By considering the above governing selection rules, we 
determined five selection-rule allowed paths for the i?i(5) two-photon transition that 
participate in the excitation process as illustrated in Fig. 2. 
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The rovibronic energies for OH A and X states calculated from the corresponding 
Hamiltonian by Luque and Crosley [7] were implemented in our simulation. The A- 
doubling and spin-orbit uncoupling effects have been incorporated into this database. 

The term energies they derived follow closely the data obtained by Dieke and Crosswhite 
[16]. A transition matrix element ( if/ v ’j\d\ tyx" r) for a transition from a lower state (v", N ") 
to an upper state (v\ N’) is given by the integral of the rovibrational wavefunctions of the 
lower and upper states together with the transition moment. A transition moment can be 
either an electronic transition moment d e {r ) (for an A-X) or a vibrational transition 
moment d v {r) (for an X-X or A- A). As for the electronic transition moment, we used the 
polynomial representation d e (r ) = «o + ayt + ayr 2 + ciyr 3 + ••• (with intemuclear distance 
r) following Ref. 17, and coefficients given therein. The vibrational transition moment is 
represented as d v (r) = er with the elementary charge e [18]. The wavefunctions (thus, the 
transition matrix elements) for all quantum states considered here were calculated by 
solving the radial Schrodinger equation using the potential energy curves for OH X and A 
states [19]. The potential curves were calculated by the first-order semiclassical 
Rydberg-Klein-Rees (RKR) method with the rotational centrifugal term [20] as shown in 
Fig. 1 . The RKR method is a widely accepted approach for the accurate determination of 
potential energy functions for diatomic molecules given the vibrational and rotational 
term energies as input parameters [21]. In the RKR calculation, the so-called Kaiser 
correction, by which higher-order semiclassical effects are taken account in some level, 
was applied to obtain further accurate potentials for the hydrides whose reduced mass is 
small [20]. A steep inner wall of the potentials was smoothed over irregular or non- 
physical behavior by the a-e' br + c function [20]. The v-dependences of G v and B v were 
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represented by the Dunham expression and its coefficients for OH X and A states [17]. 

The RKR and wavefunction calculations do not include the influence of the A-doubling 
and spin-orbit splitting (both, however, are considered in the term energies), however, 
this does not significantly impact on our spectroscopic simulation outcome as also shown 
by Luque et al [7] in their spectral simulation of OH LIF spectra. 

Finally, the energy terms and matrix elements for the two-photon transitions via the 
53 different intermediate states (0 < v < 10 for X state and 0 < v < 8 for A state) were 
accounted for in the current simulation. Note that we ignore the bounded high-lying B 2 Yf 
and C 2 S + excited states of OH and assume that they make no contribution to the two- 
photon absorption rate since their intemuclear equilibrium distances are about twice those 
of A and X [9]. 

Cross section and intermediate states of entangled-photon OH absorption. — We 
now demonstrate the unique and remarkable behavior of the entangled-photon absorption 
cross section of OH based on the theory and the molecular database described above. In 
the theoretical implementation of this technique, we follow the basic procedure described 
by Saleh et al [2]. For our simulations, we assume a monochromatic pump with the 
frequency (Op = Sf— £j (= 32,605.94 cm' 1 or 306.77 nm) and parallel-polarized entangled 
photons with A e - 10' 6 cm 2 [1,2]. The final-state linewidth function is assumed to be 
Lorentzian and the final-state radiative lifetime derived from the spontaneous emission 
coefficient is 71 1 .6 ns [7]. 

Figure 3a shows the non-monotonic variation of the entangled two-photon absorption 
cross section of OH as a function of the entanglement time T e , for the degenerate two- 
photon case. In this Figure, we also show the monotonic behavior of the classical two- 
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photon absorption cross section. In contrast with the classical absorption, a cyclic 
oscillation of the entangled-photon cross section a e is clearly visible with large variations 
that arise from the exponential term containing T e in Eq. (2). The fact that the variations 
in <j e can exceed two orders of magnitude (for T e < ~50 fs) requires a T e tuning technique 
with femto-second precision to achieve the highest level of absorption. From Fig. 3a, we 
can also see the nonlinear variation in the depth of <j e (inset right above in Fig. 3a) - this 
is because the curve is a summation of the trigonometric (exponential) functions with a 
wide range of frequency difference depending on A* . The non-degenerate absorption 

rate of OH is shown in Fig. 3b. The behavior of the non-degenerate cr e is highly 
nonmonotonic and unpredictable (without a theoretical simulation) compared with the 
degenerate case. The overall frequency of the oscillations in <j e decreases with increasing 
pump photon energy difference. It is worthy to note that a real case absorption cross 
section would likely be ensemble of degenerate and non-degenerate a e over a range of 
pump degeneracy since the SPDC spectrum is broad band, and it is not practical to use a 
narrow-band (<0.5 nm) interference filter to approximate a monochromatic pump source. 

The Fourier transform of the weighted-and-averaged entangled-photon absorption 
cross section of OH calculated from Eq. (4) is shown in Fig. 4. In order to obtain such a 
high-quality Fourier spectrum that is able to show very small peaks in Fig. 4a, 
appropriate T e and r e tunings are required over a certain range with a certain resolution. 

A reasonable range of the entanglement-time integration is found to be 20 ps for the 
current molecular system after considering the largest value of 1/| sj- £k\ (= 5.3 ps) with a 
multiplication factor of ~4X. A minimum value of the entanglement time T"" n (or 
starting value of T e tuning), which is equivalent to the maximum range of the relative- 
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path-delay r' rax , is estimated to be 67 ps when we apply the relationship that the best 
possible spectral resolution is governed by the inverse of [2]. We seek a 0.5 cm' 1 
spectral resolution for Fourier analysis (i.e. 1 IT™ ~ 0.5 cm' 1 ). The spectral range of the 

Fourier transform (Of is determined by the sampling number N s over the range of r e from 
0 to r e max , i.e. cof > nJ2t™ x . By considering a spectral range required to cover the 
possible large energy mismatches as s f —e i —co p j2 = 16,303 cm' 1 , a reasonable N s is 
found to be 2 16 , thus a r e tuning increment is r t max / N s = ~1.0 fs. 

To demonstrate the dependence of the absorption cross section Fourier spectrum on 
the T e tuning range, two spectra are calculated for the same molecular system but with 
short (3 ps) and wide (45 ps) ranges of z e at the same sampling number N s . As clearly 
shown in Fig. 5, a wider range of the relative-path-delay tuning brings a higher resolution 
for detection of peaks, thus energy levels. 

From the above parametric analysis of the entanglement tuning aspects, we found 
that a 1 0 to 1 00 ps range is required for tuning the entanglement time T e , and about 1 fs 
resolution for the tuning increment are required in order to obtain a high-resolution (0.5 
cm" 1 ) energy spectrum of the OH A—X system. 

The virtual-state spectroscopy of OH is explicitly realized in Fig. 4, which shows 
nine distinct observable peaks over three different regions of Ay with various amplitudes. 
Each peak is assigned to one of the intermediate states involved in the transition. From 
Eq. (3) for the degenerate pump case (co° k = a> p /2), the intermediate state energy is given 

by Sj — ±A j + (f, + co p / 2) . We then need a priori knowledge of the molecule’s lowest 
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intermediate state energy level [2], which is 768.36 cm’ 1 for the current system, and 
whose energy mismatch |a ; | (16078.82 cm' 1 ) is marked by an arrow in Fig. 4b. In order 
to identify the corresponding intermediate states from the peak location, we use the 
relation e j = + A ; + [s i +oj p / 2) , which should be applied to peaks whose | A y | value is 
larger than that of the lowest intermediate state. Since the amplitude of the Fourier 
transform is proportional to \A\ , only intermediate states whose transition matrix 

elements (or transition probability) Dj are relatively large appear in the spectrum. By 
reading frequencies of the peaks in Fig. 4b we are able to identify four major intermediate 
states (labeled as v ~ viii) which make a significant contribution to the absorption cross 
section. It is found that these major intermediate states are not those that lie midway 
between v" = 0 and v' = 0 (such as v" = 5 level of X 2 Hl), but they are the low-lying levels 
(v = 0) in X n, and A X. Hence, this technique is capable of identifying the exact two- 
photon paths that are dominant for the rovibronic transition. Furthermore, a remarkable 
observation is that an entanglement-induced two-photon molecular excitation from | y/vr) 
to \yvj) of OH can occur through a virtual state whose energy eigenvalue is larger than 
the incident pump energy (32,605.94 cm’ 1 ) as evident in Fig. 4c. Here, the frequency 
peak at A . = 6,068 cm’ 1 is interpreted as A 2 !, (v = 1 , N = 5), whose energy is 35,91 1 .73 

cm’ 1 . Note that this apparent violation of energy conservation can be explained because 
of the fact that this phenomenon arises from a coherent summation of the two-photon 
terms through the effects of fourth-order quantum coherence [2], which allows a virtual- 
state transition-path-scan to be performed in entanglement time- space. Specifically, if we 
examine Eqs. (4) and (5), we see that there is no limitation for the values of the energy 
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mismatch A j , while the energy conservation for the initial-to-final transition is implicitly 

ensured by the ^-function in Eq. (4). Thus, the conservation of energy is not actually 
violated. This unusual but distinct result follows a similar observation by Saleh et al [2] 
for atomic hydrogen. 

In the present non-classical technique, the virtual-states of the OH molecule in the A- 
X system are observable through a transition-path-scan which reveals individual 
intermediate-state information; such information is impossible to obtain using a classical 
two-photon spectroscopy [9]. While there are challenges to an experimental realization 
of this technique, the experimental implementation in principle is achievable if certain 
experimental parameters described above are satisfied. The new and unique features of 
entangled biphoton virtual-state spectroscopy as applied to OH demonstrate that this 
technique may be useful for spectroscopically probing difficult-to-access states of other 
molecular systems. 
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Figure Captions 

Fig. 1 The potential energy curves for the electronic ground {X) and excited (A) states of 
OH calculated using the RKR method, showing a model of a two-photon absorption 
transition s, — » q-via an intermediate (virtual) state Sj and the energy mismatch Ay 

Fig. 2 Schematic diagram of the selection-rule allowed two-photon OH rovibronic 
transitions for the R\(5)X-A absorption. The A-splitting is omitted here. 

Fig. 3 Entangled two-photon absorption cross sections for R\(5) transition in molecular 
OH A-X system over a range of the entanglement time, (a) Degenerate pump. Solid line : 
with nonclassical (entangled) light source. Dashed line : with classical light source, (b) 
Non-degenerate pump. Solid line : photon-pair frequency ratio is 1 :3. Dashed line : 
photon-pair frequency ratio is 1 :7. 

Fig. 4 Energy spectrum obtained by taking the Fourier transform of the averaged 
absorption cross section for the entangled two-photon OH A-X R\(5) transition as a 
function of the energy mismatch for an entanglement-time integration range 67 < T e < 87 
ps and a relative-path-delay tuning range 0 < r e < 67 ps with a tuning increment 5r e =1.0 
fs. The peaks are labeled by f|, f 2 , and F for Xm-m, X 2 H \ /2 , and A 2 Z state, respectively 
with (v, N) notation. 

Fig. 5 Dependence of the energy spectrum on the tuning range of r e for the same 
molecular transition in Fig. 5. Solid line : 45 <T e < 65 ps, 0 < r e < 45 ps, 8z e = 2.7 fs. 
Dashed line: 3 < T e < 23 ps, 0 < T e < 3 ps, 8r e = 0.1 8 fs. Refer Fig. 4 for the labels. 
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Fig. 3 Entangled two-photon absorption cross sections for i?,(5) transition in molecular 
OH A-X system over a range of the entanglement time, (a) Degenerate pump. Solid line'. 
with nonclassical (entangled) light source. Dashed line : with classical light source, (b) 
Non-degenerate pump. Solid line', photon-pair frequency ratio is 1:3. Dashed line : 
photon-pair frequency ratio is 1:7. 




Fig. 4 Energy spectrum obtained by taking the Fourier transform of the averaged 
absorption cross section for the entangled two-photon OH A-X R { ( 5) transition as a 
function of the energy mismatch for an entanglement-time integration range 67 < T e < 87 
ps and a relative-path-delay tuning range 0 < r, < 67 ps with a tuning increment 8z e =1.0 
fs. The peaks are labeled by f„ f 2 , and F for X 2 U V2 , X 2 Yl m , and A 2 Z state, respectively 
with (v, N) notation. 







Fig. 5 Dependence of the energy spectrum on the tuning range of z e for the same 
molecular transition in Fig. 5. Solid line : 45 <T e < 65 ps, 0 < z e < 45 ps, Sz e = 2.7 fs. 
Dashed line : 3 < T e < 23 ps, 0 < r e < 3 ps, 8z e = 0.18 fs. Refer to Fig. 4 for the labels. 




